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INTRODUCTION 



The Standard Model (SM) top quark couplings are uniquely fixed by the principle of gauge invariance, the structure 
of the quark generations, and a requirement of including the lowest dimension interaction Lagrangian \T]. Within 
the SM the top quark is considered as a point-like particle. It should be stressed, that within the SM all top-quark 
production properties and decays are evaluated with high accuracy without any phenomenological parameters. The 
total cross section production as well as the differential distributions are calculated with 0(10%) accuracy 1 1 1. The 
top quark decays through one decay channel, t bW^ (other decay channels have very small branching ratios, 
less then O(10~^)). Due to a very small life-time (rt ~ lO^^'' sec, Tt ^ 1/Aqcd) the top-quark decays long 
before it can hadronize. Therefore, it is unlikely to expect the formation of top-hadrons, T{tt)- or M{tq)-mesons 
and A(igg)-baryons. The search for anomalous (i.e. non-SM) interactions is one of the main motivations for 
studying top-quark physics. In addition, very often the events with the top quarks are backgrounds to new physics 
that we hope to discover. 

The physics of the top quark at the LHC has been studied in great detail 1 1 1, including in many cases a realistic 
simulation of the detectors (see talks at this Workshop). The goal of this presentation is to give a short review of 
Monte Carlo (MC) generators that provide a simulation of the processes with top-quark production and decays (the 
detailed considerations could be found elsewhere fT"2^). 

1 GENERAL SCHEME OF GENERATOR 

A general scheme of complete event simulation implies the evaluation of the hard process (the cross section value, 
the incoming and outgoing particle's momenta and colours), then the evolution the event through a parton show- 
ering and the hadronization of the coloured products of shower, followed by the decay of the unstable particles. 
As a result the event is described by the momenta of the final hadrons, leptons and photons and positions of 
their decay vertexes. Typically such information includes also the characteristics (momenta, colours, KF-codes, 
mother's and daughter's relations) of all intermediate partons (quarks, gluons, gauge bosons, unstable physical 
particles, etc) that allows to trace-back the history of particles production inside the event. Any such generator 
using an acceptance-rejection method (Von Neumann) returns weighted event. The most popular complete event 
generators, like HERWIG lO), PYTHIA |4|, ISAJET 0, and SHERPA include all such steps of simulation. 

The processes with top-quarks and their realization in MC generators are described in the following. The diagrams 
describing top-pair and single top production (at leading order - LO) are presented on Figs[0|2l 

'00000' - 

gg tt qq ^ tt 

Figure 1 ; tt production (QCD) 
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Figure 2: single top production (electro- weak) 

In almost all cases one needs to know the information not only about the final top-quark, but also about its decay 
products. The natural and correct way to include decays is to calculate an exact amplitude with the initial partons 
and final state particles. However, the full amplitude typically consists also of many diagrams without intermediate 
top-quarks. For examples, the full amplitude of the process gg — > tt with the both final -bosons decaying into 
electrons has many other diagrams without <-quarks. In particular, there is at least one diagram with two Z-bosons: 
gg bbZZ with one Z decaying into e+e^ and second Z — > I'ei^e- Typically, all such non-resonant diagrams 
(i.e. without top-quark) will give a very small contribution to "resonance" region, where when the invariant masses 
of the final particles are very close to i-quark and VF-boson masses (e.g. M{he^Ve) ~ mt and M{e^Ve) ~ Mw)- 
Therefore, for this case it is useful to introduce the Narrow Width Approximation (NWA). It means that one take 
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into account only diagrams with intermediate top quarks and consider all i-quarks as "on-shell" particles: 

As a result one can factorize the production and decays of the top quark. After that, typically there are three 
important issues: i) how to include top decays ?; ii) how to reproduce Breit-Wigner resonance shape ?; iii) how to 
include top spin (polarization) ? 



2 TOP DECAYS 
2.1 Decays 

Within the SM in 99.9% of the cases the top-quark decays into the bW final state: t hW, W ff. The 
matrix element describing this decay is well known fl\: 

In fact, this is the only one decay channel that is included in almost all generators. Other interesting SM channels 
(e.g t — > bW*Z* — > hf f't^l") are still not available at generator level. At the same time few other decay 
channels (due to non-SM physics) are included in generators. In particular, they are: 

o SUSY: t bH+ (almost all packages) 

o FCNC: t qg, qZ, qj (TopReX 1201 1 

o t ^ bW{-^ ff) with anomalous interactions {V + A and tensor couplings, ONETOP 1251 1 



2.2 The Breit-Wigner resonance shape 

The Narrow Width Approximation (NWA) assumes that all top quarks have the same default mass (mt). There 
are two approaches, that allow to reproduce the Breit-Wigner resonance shape due to the finite decay width of the 
top-quark. 

With the "smearing-mass" method the matrix elements are calculated for a default mt mass f4l . Then for each 
event the Breit-Wigner distribution 

/("^t) °^ 7"^1 2^2 I T^2 2 

[mt — rn^Y + T^m^ 

is used to generate separate masses (fht) for top-quarks. The momenta and energies of all final-state particles are 
re-evaluated. For example, for ti pair production with new rhi{t) and m2(t) masses one has: 

^ ± ml{t) =F 7712 (t) 

2y/I 



tt : Em = Em = ^/I/2 =^ E,,, = "'^""rr ^ W 



This method certainly is not unique, but normally provides a sensible behavior It can be used for event with any 
number of top quarks in final state (t, tt, ...). 

In another proposed method (modified NWA 1251 ) a new rht for i-quark in the event is generated by using of 
Breit-Wigner distribution. Then, the squared matrix element (|^/p) is calculated with this new mt- This method 
can be used only for the process with one top quark in the final state (i.e. not for t< production). 



2.3 How to include top polarization 

Since the top quark decays before hadronization, its spin properties are not spoiled. Therefore spin correlations in 
top production and decays is an interesting issue in top-quark physics. The simplest way to include top polarization 
employs the helicity amplitude technique |7|. An equivalent (and sometimes more simple for practical usage) 
method [SJ is realized in the TAUOLA [9J package. In this method the matrix element squared |Afp can be 
represented in the "factorized" form. 
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Mp{ab -> tX) 
Figure 3: an equivalent method (TAUOLA) 
For the process with one final top-quark (see Fig.|3j the full amplitude squared has the form: 



\M{ab-^ F + X) 



6{p^t - m^^) X \M^{ab ^ t + X)f {1 + v,h,) |A/^(i -> 



\Mp{ab -^t + X)\^ ^ |Af^(a6 ^ t + X)\^{1 + {vs)) 
\Mn{t^F){'^\MUt^F)\'{l + {hs)) 

where the amplitudes describe the production and decay of unpolarized top-quark, s is the t-quark spin 

i{pts) — 0), and v^, /i^ are the so-called "polarization" vectors L8J. 

For example, the spin-dependent matrix element squared describing top-quark decay t bt^y can be re-presented 
in the form as follows: 



|M|2 cx 



{PbPi,){ptP() 



mt(pis) 
iptPi) 



(ptPi) 



(5) 



where n| is the direction of ^+ momentum in t-quark rest frame. 



3 MATCHING ALGORITHMS 

3.1 Matching Parton Showers and Matrix Elements 

Recently a substantial progress has been achieved in the simulation of processes with multi-jet final states I10II12I . 
Indeed, the description of multi-jets obtained from the shower evolution is inaccurate, since hard radiation at large 
angle is suppressed by the angular ordering prescription. The available generators (ALPGEN, CompHEP, Madevent, 
etc) provide a generation of top production with up to 6-jets II llll3lfT4lll5l . Due to the complexity of the matrix 
elements(ME) evaluation for these many-body configurations, calculations are normally available only for LO 
cross sections. 

Multi-parton events generated using the exact leading order ME generator can be consistently evolved into multi- 
jet final states by means of a shower MC. The main problem is how to eliminate the double counting, where jets 
can arise from both the higher-order calculation and from the hard emission during the shower evolution. 

There are several approaches to this problem, aiming at different levels of accuracy. The first one ("matrix-element 
correction" technique 1 16 1) corrects the approximate ME for the emission of the hardest gluon in a given process 
by using the exact LO ME. The second is known as CKKW 1 17 |; its goal is to implement multi-jet ME corrections 
at the leading, or next-to-leading logarithmic level. This method assumes a separation of the multi-jet phase-space 
into the domains covered by the ME calculation and the domains covered by the shower evolution. Then by means 
of Sudakov re-weighting the ME's weights to reproduce the probability of an exclusive N-jet final state from the 
inclusive parton-level N-jet rate. This allows to add parton-level event samples of different jet multiplicity. 

An additional prescription is proposed in 1 1 1 1 ("MLM prescription"). In this approach the partons from the ME 
calculation are matched to the jets reconstructed after the perturbative shower Parton-level events are defined by a 
minimum Et threshold i?™™ for the partons, and a minimum separation among them, Ai?jj > Rmin- However, 
no Sudakov re-weighting is applied. Rather, events are showered, without any hard-emission veto during the 
shower. After evolution, a jet cone algorithm with cone size Rmin and minimum transverse energy E™™ is 
applied to the final state. Starting from the hardest parton, the jet which is closest to it in (ry, 0) is selected. If the 
distance between the parton and the jet centroid is smaller than Rmin, the parton and the jet match. The matched 
jet is removed from the list of jets, and matching for subsequent partons is performed. The event is fully matched 
if each parton has a matched jet. Events which do not match are rejected. For events which satisfy matching, it 
is furthermore required that no extra jet, in addition to those matching the partons, be present. Events with extra 
jets are rejected, a suppression replacing the Sudakov re-weighting used in the CKKW approach. Events obtained 
by applying this procedure to the parton level with increasing multiplicity can then be combined to obtain fully 
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inclusive samples spanning a large multiplicity range. Events with extra jets are not rejected in the case of the 
sample with highest partonic multiplicity. Fig.|4]pi"esents an illustration of above-described examples II 1II12I . 




Figure 4: Few examples of matching. The solid line corresponds to parton from a hard process, while the dashed 
one presents a parton emitted by the shower. (A) Njet = Npm-t = 3, the event is matched and kept. (B) Njet = 
Npait = 3, but A'^matchod = 2, the event should be rejected. (C) iVmatchod = A^part = 3, but Njet = 4, the event is 
matched and kept for inclusive sample, but should be rejected for exclusive samples. 



3.2 Matching with and additional b-quark 

A very simple procedure is proposed for the process of i-channel single-top production 1191 . where one has two 
complementary subprocess (see Fig.|5j: 

"2 ^ 2"qb q't and "2 -> i" qg q'tb (6) 

■t d'mm b 



t-channel 

Figure 5: Two subprocess of the single top production. 

This matching procedure is based on the transverse momentum (pt(^)) of additional 6-quark II19I . 

The matching procedure starts from simulation of the events with t-quark production due to 2 — + 2 process. Such 
an event has an additional 5-quark generated by MC showering generator (i.e. PYTHIA or HERWIG). Any event 
from this sample is accepted if the transverse momentum of the additional 6-quark is smaller than some value pq. 
The second sample consists of events from 2^3 process. The event from this sample is accepted if the additional 
6-quark from a hard process has (b) > Po- Thus, the resulting (total) sample of t-quark production events is the 
sum of two contributions: 

A^(pp ^ tX),_channel = N^^^^^pp ^ tq'; pT(b) < po) (7) 
+ N^^^^\pp^tq'b;pT{b)>Po) 
The corresponding pT(&)-distributions are shown in Fig.|5] 

This method of generation of i-channel single-top production is realized in SingleTop 1191 and TopReX 1201 codes. 



4 NLO CORRECTIONS 

Recently a substantial progress has been achieved in the calculations of the single-top production at NLO II24II25I 
1261 1271 . The radiation effects are included in the initial and final states, as well as into decays. These generators 
(like ZTOP f^Tl, ONETOP f25l and the code based on MCFM generator 1261) provide a generation of total cross 
section and differential distributions. 

One really needs to include NLO corrections in ME generators. Indeed, only shower MC's provide a representation 
of the final state complete enough to allow realistic detector simulations. On the other hand, the inclusion of NLO 
effects in fixed-order ME MC's leads to distributions which are not positive definite. The naive introduction of 
NLO matrix element could lead to double counting, since as shower MC generators already incorporate part of the 
NLO effects (a real emissions, as well as virtual effects included in the Sudakov form factors). This problem was 
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Figure 6: pt (b) distributions of the additional 5-quark production in "2 2" (the dotted curve) and "2 3" (the 
dashed curve) processes. The soHd histogram is the sum of these two contributions. The vertical line corresponds 
to parameter value — 10 GeV. 



successfully solved (see the article fW\, where it was shown how this merging can be done very effectively). This 
method is implemented in MC@NLO generator II23I describing heavy-quark pair, Higgs, DY and gauge boson pair 
production. 

The inclusion of NLO corrections in the shower MC guarantees that total cross-sections generated by the MC 
reproduce those of the NLO ME calculation, thereby properly including the K factors and reducing the system- 
atic uncertainties induced by renormalization and factorization scale variations. At the same time, however, the 
presence of the higher-order corrections generated by the shower will improve the description of the NLO dis- 
tributions, leading to departures from the parton-level NLO result. This is shown, for example, in Fig. which 
shows the pr spectrum of a tt pair resulting from the pure NLO calculation, from the LO shower, and from the 
MC@NLO improvement. At large px, a region dominated by the NLO effects, MC@NLO faithfully reproduces the 
hard, large-angle emission distribution given by the NLO matrix elements. At small pT, a region dominated by 
multiple radiation and higher-order effects, the MC@NLO departs significantly from the NLO result, while properly 
incorporating the Sudakov re-summation effects only available via the IS shower evolution. 




logio(PT"/GeV) 



Figure 7: Transverse momentum distribution of top quark pairs using three different approaches: the LO HERWIG 
MC, the parton-level NLO MC, and the merging of the two into MC@NLO. Figure from 1221 . 



Recently, NLO corrections to the single-top production processes [t- and s-channels) are included at MC@NLO 
event generator 1271 . The comparison MC@NLO and HERWIG results for the single top production are given on 
Fig.|8l Note, that the both generators give very close results for the highest pT jet distribution (the left histogram 
on Fig.|8j, while MC@NLO predicts much harder pT-spectrum for fe-jet (not from top decay). 



5 EVENT GENERATORS 

This section presents a brief description of the generators. Only few basic comments and the list of included 
top-quark production process are given for each entry. 
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Figure 8: The transverse momentum distributions for the highest px jet and 6-jet (not from top decay). 



5.1 Full event simulation packages 

These packages provide a full event simulation including the hard process generation, showering and hadronization 
with subsequent decays of the unstable hadrons. 

•herwigEI 

contains a wide range of Standard Model, Higgs and supersymmetric processes. It uses the parton-shower approach 
for initial- and final-state QCD radiation, including colour coherence effects and azimuthal correlations both within 
and between the jets. HERWIG is particularly sophisticated in its treatment of the subsequent decay of unstable 
resonances, including full spin correlations for most processes. 
Processes included: tt, single top (t, s channels), ttH, Ztt, gb tH+ 



• PYTHIA |4| 

is a general-purpose generator for hadronic events in pp, e+e^ and ep colliders. It contains a subprocess library 
and generation machinery, initial- and final-state parton showers, underlying event, hadronization and decays, and 
analysis tools. 

Processes included: tt, single-top (t, s channels), ttH, gb tH^ , no spin correlations 

• ISAJET 13 

is a general-purpose generator for hadronic events. ISAJET is based on perturbative QCD plus phenomenological 
models for parton and beam jet fragmentation. 
Processes included: tt, no spin correlations 

•SHERPA |6l 

is a new multi purpose event generator with a powerful matrix element generator AMEGIC++ 



5.2 Tree level matrix element generators 

Such packages generate the hard processes kinematic quantities, such as masses and momenta, the spin, the colour 
connection, and the flavour of initial- and final-state partons. Then such information is stored in the "Les Houches" 
format 1291 and is passed to full event simulation generator (like PYTHIA or HERWIG). 

•ALPGEN JTH 

is designed for the generation of the SM processes in hadronic collisions, with emphasis on final states with large 
jet multiplicities. It is based on the exact LO evaluation of partonic matrix elements, as well as top and gauge 
boson decays with helicity correlations. The code generates events in both a weighted and unweighted mode. 

Processes included: ti + up to 6jets, single top: tq, tb, tW, tbW (no extra jets), tttt + up to 4jets, ttbb + up to 
4jets, Htt + up to 4jets, W/Ztt + up to 4jets 

•GompHEP IT3l 

CompHEP computes squared Feynman diagrams symboUcally and then numerically calculates cross sections and 
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distributions. The output event's information in the "Les Houches" format [29 1. 
Processes included: ti, single top, (? tttt, ttbb, ) W/Ztt, spin correlations are included 

•MadEvent [WJ 

is a multi-purpose, tree-level event generator, which is powered by the matrix element generator MadGraph 1151 . 
MadGraph automatically generates the amplitudes for all the relevant subprocesses and produces the mappings for 
the integration over the phase space. 

Processes included: tt + up to 3jets, single top, ttbb + up to Ijet, Htt up to 2jets 

• MC@NLO [23] combines a Monte Carlo event generator with exact NLO calculations of rates for QCD pro- 
cesses at hadron colliders. 

Processes included: tt, single top (t- and s-channel) 
•AcerMC |21 1 (see talk of B.Krersevan in this conference). 

AcerMC is dedicated for generation of the Standard Model background processes in pp collisions at the LHC 
Processes included: tt, single top, tttt, ttbb, W/Ztt, spin correlations are included 

•SingleTop 1191 generator is based on the CompHEP 1131 package. 

processes: t-channel single top production (2 ^ 2 h- 2 3, spin correlations are included 
•TopReX|20| 

provides a simulation of several important processes in hadronic collision, not implemented in PYTHIA. Several 
top-quark decays channels are included: the SM channel {t qW'^, q — d,s, b), 6-quark and charged Higgs 
{t bH~^) and the channels with flavor changing neutral current (FCNC): t — + u{c) V, V — g, 7, Z. The 
implemented matrix elements take into account spin polarizations of the top quark. 

Processes included: gg (qq) tt, single-top production {t-, s-, and iVF-channel), qq' H^* tb, qq — + 
W*/Z*QQ with W*/Z* -> // and Q = c,b, t, gu{c) t bW (due to FCNC). 

• MCFM 1281 included the matrix elements at next-to-leading order and incorporate full spin correlations. 
Processes included: tt, single top (t- and s-channel), Htt, W/Ztt 

•ZTOP 1241 code includes the full NLO-corrections to single top production (t- and s-channel). 

•ONETOP 1251 code include full NLO-corrections to single top production [t- and s-channel) and top-quark 
decay. 

CONCLUSIONS 

Recently substantial progress has been achieved in the implementation of new ideas concerning top events simula- 
tion and in developing event generators for top quark production and decays. 

In particular, several new processes describing top-quark production and decays are included in the generators 
providing complete event simulation. Many event generators provide generation of top-quark processes with spin 
correlations. Tree level generators make possible the generation of top production processes with additional multi- 
jets in the final states. The codes with full NLO corrections to single-top production processes are available. A few 
generators include also processes with top-quark production and decays due to interactions beyond SM. 
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